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Abstract 

The motion of an elastic body over a rough surface is simulated in one 
dimension using Newtonian dynamics. We drive the body by prescribing 
a constant velocity. We extract the friction force from the work done by 
the applied force in order to maintain the motion of the body. A positive 
amount of work is needed because energy is converted into the internal 
elastic degrees of freedom of the body. The velocity dependence of the 
friction force is studied as function of the roughness of the surface and the 
elastic properties of the mobile body. We pay special attention to the low 
velocity limit. It is shown that the friction force goes to a finite value as 
the velocity goes to zero when the substrate is sufficiently rough to induce 
elastic instabilities. For smooth substrates, where elasitc instabilities are 
absent, the friction force goes to zero when the velocity goes to zero. 

PACS numbers: 46.30.Pa, 81.40.Pq, 68.10.Et, 05.70.L 

1 Introduction 

The present paper represents an attempt to understand how the the velocity 
dependence of the friction force is related to the roughness of the substrate and 
to the elastic properties of the sliding material. 
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The role friction plays in our everyday lives, and many technological appli- 
cations, is immense and yet there still does not exist a successful quantitative 
theory that explains this phenomena. The variation of the frictional force with 
the sliding velocity is interesting in a practical engineering sense where the aim is 
to control or eliminate friction, but it is also of relevance to workers in a number 
of other fields. Experimentally, dynamic friction is important in the study of 
stick-slip motion [|TJ, for example. Theoretically, friction is relevant, for instance, 
to models of earthquake faults @, |^, || and railway wheel squeal 0. In such prob- 
lems a specific form of the velocity dependence of the friction force is postulated 
on phenomenological grounds. The friction law may be nonlocal and nonlinear 
[D, [?]], local and nonlinear p], [| ^| or simply local and linear ||. 

Cieplak, Smith and Robbins[[| called the linear relationship // r (f) oc v be- 
tween velocity v and the friction force ff r for viscous friction in their beautiful 
study of the microscopic origin of the friction force. If lim„^ ffr 7^ they refer 
to the behaviour as static friction. They found a change from viscous to static 
friction as the corrugation of the substrate potential is increased. We believe that 
this change in velocity dependence is related to the onset of elastic instabilities 



similar to those described in Ref. [|I(], [Ll| |12| [13], |14|| . When the local curvature 
of the substrate potential becomes sufficiently large the static force balance equa- 
tion for the individual atoms supports multiple solutions. The abrupt jump from 
one of these solutions to the other leads to a finite dissipation even in the limit 
v — > 0. Cieplak et al. discuss at the end of their paper Ref. || the reason why 
static friction is more often encountered than static friction. We suggest that 
the reason migh be that a slight roughness at the interfaces will produce elastic 
instabilities and therefore a finite friciton force in the limit of vanishing velocity. 

Shinjo and Hirano have observed friction in a Frenkel-Kontorova type of model 



with added kinetic energy terms [jlgl . The friction force they observe in the limit 
of low velocity is non-zero when the coupling in the system is greater than a 
critical value, this being the Aubry transition point jl6|. Shinjo and Hirano also 
found a vanishing friction force in the limit of high velocity. Their findings imply 
that elasic instabilities are important for the existence of a finite static friction 
force. This is in agreement with the work presented in the present paper. Sokoloff 
has also made the point that local jumping motion is essential for the occurrence 



of a dynamical threshold force in a microscopic model of idealized surfaces 
He suggests that this is probably also applicable to macroscopic surfaces if we 
think of these surfaces as consisting of two elastic media which interact more 
strongly at some points on their surfaces than others. Our macroscopic model 
shows behaviour which is consistent with this idea. 

Sokoloff has performed simulations which give clear evidence for a transistion 



from friction free motion to a dissipative regime as function of system size [18] 
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The transition can be understood in terms of Chirikov's overlap criterion [I8[, |T9[ . 
This criterion states that the frequency of the driving force (for instance the 
wash board frequency of the substrate potential) must overlap with one of the 
eigenfrequencies of the elastic modes in order for energy to be dissipated into the 
phonon system. Since the spacing between the eigenfrequencies decreases with 
increasing system size the overlap criterion is likely to be fulfilled as the size of the 
system is increased. This is what Sokoloff has demonstrated in fact happens. The 
elastic instability induced friction mentioned above can be seen as an alternative 
mechanism for inducing dissipation. Chirikov- Sokoloff resonance mechanism is 
of dynamical origin, whereas the elastic instability mechanism originates in the 
possible multivalued character of the static force balance. 

A wealth of experimental work has led to the qualitative aspects of friction 



being well defined (for a review see [2(| fJT], j22|). Recently, interest has 



concentrated on a microscopic description and understanding of the origins the 



friction force O, U5L 23, Experimental techniques such as atomic force 



microscopy (AFM) ]25], ^] and friction force microscopy (FFM) allow the 



processes of friction to be studied in three dimensions and ah initio molecular 
dynamics simulations have been used to study atomic scale mechanisms of energy 
dissipation ||. Work by Sokoloff [[| has found the velocity dependence of the 
friction force in a model where the internal degrees of freedom of the moving 
system are subject to linear damping described by a damping coefficinet 7. The 
obtained velocity dependence of the friction force remains unchanged in the limit 
of 7 -»• 0. 

We have previously studied the limit v — > 0. In this limit we were able to 
show that in order to calculate the restoring static friction force perturbatively one 
must expand at least to third order in the random substrate potential WM. This 
is consistent with the observation that a finite static friction force is produced by 
the existence of non-linear elastic (or plastic) instabilities \T2\. We have studied 



the effect as well as the statistical properties of these instabilities [T3|, |28| . In these 
studies we assumed that the time evolution was controled by a set of over-damped 
dynamical equations. We did not study the mechanisms behind the damping. In 
the present paper we attempt to treat the damping in a more realistic maner. 
We study an elastic system consisting of a large number of internal degrees of 
freedom (light particles) which do not interect with the external potential. These 
degrees of freedom act the heat bath which is able to absorb the energy released 
during the instabilities. 

In the simplest phenomenological description one assumes some specific func- 
tional form for velocity dependence of the friction force. There are two main 
classes of behaviour. These are namely velocity strengthening, where the friction 
increases with velocity, and velocity weakening, where the opposite occurs, often 
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in a nonlinear fashion. A simple example of a velocity strengthening friction 
function is the traditional proportionality assumption f/ r = — rjv. However, of- 
ten one observes velocity weakening of the friction force. Velocity weakening and 
strengthening behaviour have been observed in experiments with cast iron pO 



and recently with rubber against glass |29| and serpentinite |3(| . One will expect 
strengthening to cross over to weakening for large velocities when the interaction 
time between the individual extremal point of the two sliding surfaces goes to 
zero. This is, indeed, what we observe in our simulations. 



2 One Dimensional Model 



The specific model we consider is similar to the Burridge-Knopoff spring-block 
model|T[[]. This model has been studied recently by several authors[TJ, £| ffl. Our 
study differs from these investigations in the following important way. The work 
in Ref. || |||, f§] assume a specific form for the velocity dependent friction force 
acting on the individual blocks. In the present study we use purely Newtonian dy- 
namics and study the velocity dependence of the effective friction force produced 
by dissipation of kinetic energy of the blocks into internal energy of a mechanical 
heat bath attached to the blocks. 

We have the relative motion of two macroscopic bodies in mind. We imagine 
the two interfaces to be rough on an atomic scale and focus on the interaction 
of the asperities reaching out from the two surfaces. See insert to Fig. 1. For 
simplicity we model the lower material by a static potential. The elastic de- 
formations induced in the sliding upper material are captured by the following 
somewhat oversimplified Hamiltonian. We focus our attention at only two as- 
perities which we represent by two heavy particles of mass M and position qi 
(See Fig. 1). I.e. the asperities are modeled as stiff bodies. They interact with 
a substrate potential U(qi). The induced elastic deformation is represented by 
the elastic distortion of the bar of eigenfrequency Q coupling the two particles 
together. The heat bath consituted by the phonons in the bulk of the material 
is represented by an elastic chain of N light particles of mass m and position 
x n coupled to the two heavy particles. The overall stifness of the material is 
modelled by the spring of stiffness K = MQ 2 . The Hamiltonian of the system is 
given by 

H = ^Mqf + ^Mq 2 2 + ^Mn 2 ( qi ~q 2 -L) 2 (1) 

1 1 

+-mu 2 (q 1 -xi- a) 2 + -muj 2 (x N - q 2 - a) 2 + U(q 1 ) + U(q 2 ) - q\F - q 2 F 
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N 1 1 

+ X, \^ mi: n + ^mu 2 (x n - x n+1 - a) 2 } 

n=l 1 1 



Where L = (N + l)a and u is the eigenfrequency of the springs connected to the 
light particles. We do not include any friction term in the dynamical equation of 
the model. The dynamics of the model are directly given by Newton's equation 
md 2 x/dt 2 = —5H/5x where x = q%, g 2 , • • • , and calculated numerically 



using the leapfrog algorithm [ 32] . The potential U representing the rough surface 
is produced by randomly positioning Gaussian pinning centres with density n p 
along the x-axis. Each individual pinning well U p has the form 

U p (r) = -A p eM-(r/R P ) 2 ). (2) 



We drive the model by prescribing a centre of mass velocity, v com which is 
then kept constant during the simulation. This is obtained by adjusting the 
applied force F in each time step such that the total force on the system F tot = 
2F — dU/ dq\ — dU / dq 2 is always kept equal to zero. The friction froce is identified 
as the time average of the force F . The average is taken over a time window 
during which the rate of change of the internal kinetic energy of the heath bath 
is constant. This time interval begins after a short transient period after which 
the kinetic energy of the bath increases linearly with time. We make sure to stop 
the measurement before the process becomes non-stationary due to heating of 
the heat-bath degrees of freedom. 

We used the following set of units. The mass M = 1, the length L = 1, 
and time 1/0=1. A series of constant velocity simulations were carried out 
using a range of values for the pin amplitude A p , small spring constant mio 2 and 
pin range R p . All simulations were done with a light particle mass m = 0.1, 
number of small particle = 999, and a density of pinning wells n p = 1000. The 
time step of the discrete numerical integration of the equation of motion were 
A = 0.001. All results reported were found after averaging over a number of 
realizations of the background potential. 



3 Simulation Results 



We are particularly interested in the effect of elastic instabilities. Let us briefly 
discuss the nature of these instabilites. Consider a particle of position x elastically 
coupled by a spring of stiffness A; to a moving coordinate X(t) = Vt, where V is 
a constant velocity. I.e. the elastic energy of the particle is E el = ^k(x — Vt) 2 . 
Let the particle move through a pinning well of the form given in Eq. |2|. It is 
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easy to show (See the third paper of Ref. |12|) that the force balance equation 
for the particle in the limit V — > will have multiple solutions (See Fig. 2) when 

The particle undergoes discontinous jumps (in the limit where X(t) moves for- 
ward quasistatically) when the solution of the stability equation becomes multi- 
valued. The two heavy particles of our model will influence each other during 
the instabilities. The condition for instabilities in Eq. ||| is therefore only approx- 
imately applicable to our model. 

We present in Fig. 3 simulations of the velocity dependence of the friction 
force for parameters around the threshold value given in Eq. |3|. Fig. 3a presents 
a set of data for different values of the range of the individual wells and Fig. 3b 
show the behaviour for different amplitudes of the pinning wells. 

One observe in Fig. 3a a crossover in the behaviour at low velocities for R p 
about equal to the threshold value R* = (4A p /e 3//2 fc) 1//2 . Similarly in Fig. 3b a 

change in the behaviour occur for A p around A* = ^-kR? v . In the case where 
Eq. |3] predicts the absence of elastic instabilities in the limit v — > one see that 
the the friction force is an increasing function of the velocity for small velocities. 
The intercept with the y-axis vanishes gradualy as the instabilities cease to ex- 
ist. This behaviour can clearly be identified despite the difficulties of producing 
high quality statistics in the limit of vanishing velocity. (Remember that due to 
the randomness in the position of the pining centres the local curvature can be 
stronger than the maximum curvature of a single well.) 

That a non-zero value of lim^o ff r (v) is linked to the existence of elastic 
instabilities is more clearly demonstrated in Fig. 4. In Fig. 4a we show the time 
dependence of the kinetic energy of the bath as one of the heavy particles moves 
through a single potential well. Fig. 4b shows how the change in the energy 
of the bath depends on the velocity of the system. It is unfortunately difficult 
to simulate the very low velocity limit with high accuracy. However, one sees 
that AE^in (the change in the energy of the bath) can be extrapolated to zero as 
v —> when A p < A*. 



4 Discussion 



The above observed behaviour of the friction force can qualitatively be understood 
in terms of the following simple model. Consider a single particle of position x. 
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Let the particle of mass m be coupled harmonically, through a spring of stiffness 
mwo, to a frame which is moving with velocity v. Assume furthermore that the 
particle is in contact with a heat bath which enables a damping propotional to 
the velocity of the particle to take place. The interaction with a rough substrate 
potential is represented by the force F. In the frame movng with velocity v the 
equation of motion is 

d 2 x dx 
m ~dp = ~ muJ oX - 7— + f(x(t)). (4) 

We can solve this equation immediately if we replace the position dependent 
force f(x(t)) by a time dependent force F(t). We shall chose the the following 
mean field-like form F(t) = f(vt). This is a reasonable approximation for smooth 
over-damped motion. The time averaged work-per-time involved in the motion 
is given by 

1 f T / 2 „, ,dx 



The work is readily expressed in terms of the Fourier transform of the Green's 
function _ 

rod 

X (uj) = / dte^xit) = [m{ul - u 2 ) - i^]' 1 (6) 



of Eq. |j. One obtains 

/OO fj, ) 
— (F(u)F(-u))(-iw) X (u). (7) 
-oo Z7T 

We introduced the Fourier transform of the force auto-correlation function 

poo 

(F(u)F(-u)) = / dte^(F(t)F(t + t )) t(y (8) 



The time dependence of F(t) is taken to be given by a linear super-position 
of the individual forces exerted by the randomly superpositioned pinning wells. 
A pinning well positioned at Xi is encountered at time t — Xi/v. We have 

oo 

F(t)= E fp(t~ti)- (9) 

i=— oo 

The auto-correlation of F(t) is according to Campbell's theorem given by 

POO 

(F(t)F(t + *„))*, = n p v / dt f p (t )f(to + t) (10) 



We will consider two types of time evolution of f p (t). Firstly, in the absence of 
instabilities f p (t) will be a smooth function of t. We chose to represent f p (t) by 
the continuous function 

W) - -^- ( *' 2 . (id 

-ftp 
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We have introduced the time scale A = R p /v. The second type of time depen- 
dence involves to rapid changes in f p (t) encountered when an elastic instability 
takes place. See Fig. 2b. We will represent these abrupt jumps in f p (t) by the 
following function 

f P (t) = m-a±-Q(A-t) (12) 



The force correlator is easily calcuated in both cases. When no instabilities 
exist one finds 

(F c (uj)F c (-uj)) = ™p(Vg 2 a; 2e i (A< (13) 



v- 



When elastic instabilities take place the correlator is slightly more complicated 
(F(u)F(-u)) = (F c {uj)F c {-uj)) + 8A p n p a Sm ^^ f°° dzsm{Aujz)e- z ( 14) 

UJ J— oo 



^^[^(cos(2A^) + l) 
A uj z 

2 4A 
■— (cos(2Acj) - 1) -sin(2Ao;) . 

UJ UJ 6 



It is now easy to calculate the friction force ff r = W/v from Eq. |7|. The 
friction force vanishes in both as an inverse power of v in the high velocity limit. 
In the low velocity limit one finds in the case of no instabilities 

3 [W in v A] 
fr ~ 4 V 2m 2 ^ 4 



ffr - Ixh-^JV. (15) 



I.e. the friction force goes to zero linearly with v. The presence of instabilities 
produces a non-zero value of friction force in the limit v = 

f fr (v = 0) = ^. (16) 



muj, 



2 ' 



The instabilities also change the coefficient to the term in ff r (v) linear in v. One 
note that the zero- velocity limit in Eq. |TB| does not depend on the damping co- 
efficient 7. This is to be expected since in the lmit of quasi-static motion any 
non-zero value of 7 will suffice to extract the energy relased during the instabil- 
ities. It is important to remember that the coefficient a describing the jump in 
the substrate force depends implicitly upon the stiffness of the elastic coupling 
of the moving particle. The softer the elastic coupling the larger the jumps (See 
paper 3 of Ref. [fl2]0. This does not change the qualitative significances of Eq. 
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5 Conclusion 



We have studied a one dimensional model of an elastic body moving over a rough 
surface. We focus our attention on the behaviour of the macroscopic asperities 
poking out into the rough surface. The body is moved at constant velocity. We 
find that the velocity dependent friction force always decreases to zero in the 
limit of large velocity. The value of the friction force as the velocity tends to 
zero depends on the corrugation of the substrate potential. If the curvature of 
the potential is large enough to induce elastic instabilities during the quasistatic 
motion the friciton force goes to a non-zero limit. This behaviour was called static 
friction by Cieplak et al.0. For smooth substrate potentials elastic instabilities 
may be absent and then the friciton force goes to zero with vanishing velocity. 
Cieplak et al. called this velocity dependence for viscous friction. 
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Figure Captions 



Figure 1. 

Schematic diagram of the one dimensional model used in the simulations. We 
model the relative motion of two macroscopic bodies (inset) and focus on the 
interaction between two asperities protruding from the upper body and a static 
potential, which represents the lower body. Elastic deformation is represented by 
the long elastic bar and a heat bath by an elastic chain of light particles. 

Figure 2. 

(a) The spatial variation of the force experienced by an harmonic oscillator as it 
is pulled through a single pinning centre. The slope the dashed lines is the spring 
constant of the oscillator. When the critereon [3] is met the force balance equation 
supports mutiple solutions and the oscillator will experience jumps in position, 
represented by the arrows. This gives rise to abrupt jumps in the temporal 
varition of the force experienced by the oscillator, as shown in (b). 

Figure 3. 

Simulation results showing the velocity dependence of the friction force for dif- 
ferent values of (a) the range R p of individual pinning centres; (b) the amplitude 
A p . R p * and A p * are the threshold values for the pin range and amplitude re- 
spectively. 

Figure 4. 

(a) Time dependence of the kinetic energy of the heat bath as one of the heavy 
particles moves through a single potential well. Values were taken for different 
values of the pin amplitude A p . (b) Velocity dependence of the change in the 
kinetic energy of the heat bath during the motion of one heavy particle through 
a single potential well. Values were taken for different values of the pin amplitude 
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Figure 1. Lindop and Jensen. 




Figure 2. Lindop and Jensen. 
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Figure 3(a). Lindop and Jensen 
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Figure 3(b). Lindop and Jensen 
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Figure 4(a). Lindop and Jensen 
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Figure 4(b). Lindop and Jensen. 




